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Abstract 

Using continuation methods from the integrable Ablowitz-Ladik lattice, we have 
studied the structure of numerically exact mobile discrete breathers in the standard 
Discrete Nonlinear Schrodinger equation. We show that, away from that integrable 
limit, the mobile pulse is dressed by a background of resonant plane waves with 
wavevectors given by a certain selection rule. This background is seen to be essential 
for supporting mobile localization in the absence of integrability. We show how 
the variations of the localized pulse energy during its motion are balanced by the 
interaction with this background, allowing the localization mobility along the lattice. 
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1 Introduction 



The phenomenon of intrinsic localization (collapse to self-localized states) due 
to nonlinearity in discrete systems governed by Schrodinger equations is of 
fundamental interest in Nonlinear Physics [1,2], and is the subject of current 
active experimental research in several areas like nonlinear optics [3], Bose- 
Einstein condensate arrays [4,5,6], polaronic effects in biomolecular processes, 
and local (stretching) modes in molecules and molecular crystals (see [1,7,8] 
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and references therein). Discrete Non-Linear Schrodinger equations (NLS lat- 
tices for short) provide the theoretical description of these systems, where 
pulse-like (self-localized) states are observed. 

The standard Discrete Nonlinear Schrodinger (DNLS) equation is the (sim- 
plest) discretization of the one- dimensional continuous Schrodinger equation 
with cubic nonlinearity in the interaction term, i.e., 

i$ n = -($ n+ i + $ n _i) -7|$ n | 2 $ n , (1) 

where & n (t) is a complex function of time. The first term on the right takes 
account of the dispersion and the second of the nonlinearity, the parameter 7 
is the ratio between them. For the Bose-Einstein condensate lattices dealt with 
in [4,5,6] one can think of $ n as the boson condensate wavefunction in the n-th 
(optical) potential well, and 7 would thus be related to the so-called s-wave 
scattering length [9]. The self-focussing effect of local nonlinearity balanced by 
the opposite effect of the dispersive coupling makes possible the existence of 
localized boson states in the Schrodinger representation of the condensate lat- 
tice (Gross-Pitaevskii equation). In a localized state (discrete breather) of the 
boson lattice the profile of |$ n | 2 decays exponentially away from the localiza- 
tion center. These solutions have an internal frequency, $ n = |$„| exp(iu^t), 
so that the discreteness is essential to avoid resonances with the phonon band 
and keeping localized the energy. Pinned (immobile) localized solutions of eq. 
(1) have been rigorously characterized [10] and extensively studied by highly 
accurate numerical [11] and analytical approximations. However, for exact 
mobile discrete breathers no rigorous formal proof of existence in standard 
DNLS is available nowadays although lot of works have studied these kind of 
solutions (see e.g. [12,13,14,15]). 

The translational motion of discrete breathers introduces a new time scale 
(the inverse velocity) into play, so generically a moving breather should excite 
resonances with the plane wave band expectra. In a hamiltonian system, these 
radiative losses would tend to delocalize energy and some compensating mech- 
anism is needed in order to sustain exact stationary states of breather trans- 
lational motion. To address the problem we use unbiased (i.e. not based on 
ansatze on the expected functional form of the exact solution) and precise nu- 
merical methods which allow observations of numerically exact non-integrable 
mobility, paving the way to further physical (and mathematical) insights. 

In this letter, after explaining in section 2 the basis of the numerical method 
(fixed point continuation from the integrable Ablowitz-Ladik limit [16]) and 
its relevant technical details briefly, we will discuss the structure of the discrete 
NLS breathers in 3. They are found to be the exact superposition of a travelling 
exponentially localized oscillation (the core), and an extended "background" 
built up of finite amplitude plane waves Aexp[i(kn — out)]. These resonant 
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plane waves fit well simple (thermodynamic limit) predictions based on dis- 
crete symmetry requirements. Finally in section 4 we show how the resonant 
background is seen to be an indispensable part of the solution. In this regard we 
present the mechanism through which the interaction core-background com- 
pensates the variations of the core energy (no longer an invariant of motion 
away from the integrable limit), during the translational motion. 



2 Salerno Model and Continuation Method 



The method used here makes use of the following NLS lattice, originally in- 
troduced by Salerno [17], 



i$ n = -($ n +i + $ n -i, 



2z/$J$ r 



(2) 



This lattice, though non-integrable for v ^ 0, provides a Hamiltonian interpo- 
lation between the standard DNLS equation (1), for \i = and v = 7/2, and 
the integrable Ablowitz-Ladik lattice [16], A-L for short, when \i = 7/2 and 
v = 0. The A-L model is a remarkable integrable lattice possessing a family 
of exact moving breather solutions: 



■ sinh (3 sech[/5(n — x (t))] exp[i(a(n — x (t)) + fi(t))], 



(3) 



the two parameters Uf, and Vf, are the breather frequency and velocity 



u h = Cl{t) = 2 cosh p cos a + avi, , = x (t) = — sinh f3 sin a (4) 



where — n < a < n and < j3 < 00. The equation (2) has the following 
conserved quantities, namely the Hamiltonian H and the norm iV: 



H = - A+i + $ n $ n+1 ) - 2- ]T |$„| 2 + 2— £ ln(l + /i|$„| 2 ) (5) 

n ^ n A 4 n 

7V = I^ln(l + / ,|$ n | 2 ) (6) 

A 1 n 

where $ n denotes the complex conjugate of <£>„. In what follows we will fix the 
value 7 = 2 in eq. (1) and /i + v — 1 in eq. (2), as usual. 

Perturbative inverse scattering transform [18], as well as collective coordinate 
methods [19,2,20], have been used to study moving breathers of the Salerno 
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equation (2) near the integrable A-L limit (is ~ 0). The numerical procedure 
that we explain below has the advantage of being unbiased and not restricted 
to small values of the non-integrability parameter is, at the expense of restrict- 
ing attention to those solutions (3) which are resonant, meaning that the two 
breather time scales are commensurate 2nv b /uj b = p/q (rational time scales 
ratio). A resonant (p/q) moving breather $ n (t) is numerically represented as 
a fixed point of the map M = L p T q , where L is the lattice translation oper- 
ator L({<3> n (t)}) = {$ n+1 (t)}, and T is the T b -evolution map (T b = 2it /uj b ), 
T({$ n (t)}) = {$ n (t + T b )h explicitly 

$„(t) = $ n+p (t + qT b ) for all n (7) 

Let us briefly present the numerical method. The implicit function theorem 
[21] ensures a unique continuation of a fixed point solution of M for parame- 
ter (is) variations, provided the Jacobian matrix J = D(M — I) is invertible: 
with this proviso the Newton method [22] is an efficient numerical algorithm 
to find the uniquely continued fixed point. In other words, continuation from 
a resonant A-L breather along the Salerno model is possible if one restricts 
the Jacobian matrix J to the subspace orthogonal to its center (null) sub- 
space. The center subspace turns out to be spanned by two continuous sym- 
metries of the Salerno model, namely, time translation and gauge (uniform 
phase rotation) invariances. Using Singular Value Decomposition (SVD) tech- 
niques [23] , one then obtains numerically accurate continued resonant moving 
breathers along the Salerno model until conditions for continuation cease to 
hold. A (SVD)-regularized Newton algorithm was already used by Cretegny 
and Aubry in [12] to refine moving breathers of Klein-Gordon lattices with 
Morse potentials obtained by other means. From the methodological side 
what is novel here is the systematic use of it in order to obtain the family 
of moving Schrodinger breathers of the NLS lattice (2), for different values of 
2nv b /u) b = 0, 1/2, 3/4, 1,... and a fine grid of frequency values u b and the 
nonintegrability parameter v. 



3 Mobile Discrete Breathers 

Let first start with a few remarks on immobile breathers (p — 0). Some non- 
integrable issues, that affect mobile solutions, can be shown continuing the 
immobile ones along the Salerno Model (is — 0,. . . , 1). First we remark that 
the uniquely continued solution of standard DNLS (is = 1) is equal to the 
pinned discrete breather uniquely continued from the anticontinuous limit 
(7 — > 00) [10]. Second, only inmobile breathers which are centered either at 
a site (n) or at a bond (n ± 1/2) persist; this is due to the emergence of 
Peierls-Nabarro barriers away from integrability (is ^ 0), a well-known result 
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of collective variable theory [2,20]. The breather centered at a site is stable 
while the one centered at a bond is unstable [2,24]; its energy difference is 
the Peierls-Nabarro barrier. This energy difference acts as a barrier to mobile 
breathers for travelling along the lattice; the numerical computations of this 
barrier nicely fit with collective variable predictions. 

Our main interest, however, focusses on mobile solutions, i.e. p 7^ 0. How 
are Peierls-Nabarro barriers to mobility overcome by the fixed point solution? 
Our results show clearly that the uniquely continued p/g-resonant fixed point 
for v 7^ is spatially asymptotic to an extended background, whose amplitude 
increases from zero (at v — 0) with increasing non-integrability v, superposed 
to the moving (A-L)-like core, see figure (1). In order to reveal the structure 
of this extended background, we have to pay attention to spatially extended 
solutions of the Salerno model. 



(a) (b) 




n 

Fig. 1. Instantaneous profile of a 1/1 breather with uj^ = 4.909 and Vb = 0.7813. 
Real part (a), imaginary part (b), modulus (c) and phase (d). The nonintegrability 
parameter in (2) is v = 0.05. 

The Salerno equation (2) admits extended solutions of the plane wave form, 
$ n (i) = Aexp[i(kn — cut)], provided the following (nonlinear) dispersion rela- 
tion holds: 

u = -2[1 + (1- u)\A\ 2 ]cosk-2u\A\ 2 (8) 
A p/g-resonant plane wave satisfies $„(t) = <& n+p {t + qTj,), and therefore, for 
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a p/g-resonant plane wave the following condition also holds: 




) 



(9) 



m being any integer. Equations (8) and (9) can be solved for k and one obtains 
a finite number of branches of p/g-resonant wavenumbers in the first 

Brillouin zone, — tt < kj < it. The simplest case of a unique branch for fixed v 
(as well as uj^ and p/q) and A small, is represented in a) and b) of figure (2). 
For example, for A small and u>b > 4, for any value of < v < 1 and A, there 
is a unique 1/1-resonant wavenumber branch ko(v, A). 

For the general situation where several branches kj (j — 0, . . ., s — 1) of 
resonant plane waves solve (8) and (9), the power spectrum of a background 
site n, S(u) = | Pf^ 9?(<E> n (i)) exp[iu;t\dt | 2 , reveals s peaks at the values uj 
corresponding to the resonant wavevector branches. The background is, up to 
numerical accuracy, a linear superposition of pj g-resonant plane waves, namely 



The amplitudes Aj differ typically orders of magnitude, i.e. \A \ ^> \Ai\ ^> 
\A 2 \. . ., so that only a few frequencies are dominant, for most practical pur- 
poses. One would speak of localization in /c-space to describe the extended 
background of the p/g-resonant fixed point. Once the values of uj^, t>6 and v 
are given, the "selection rule" provided by equations (8) and (9), does not de- 
termine directly the resonant wavenumbers kj, but only branches kj(A). This 
reflects the inherent nonlinearity of the NLS lattice, wherefrom the frequency 
of the plane wave depends on both wavenumber and amplitude in equation 
(8). Along the parametric continuation path the fixed point "adjusts" the 
planewave content (kj) of the background, so that it remains p/g-resonant 
under the changes in the amplitudes of the background plane waves (Aj). 



4 Background relevance to Mobility 

Along the Newton continuation path to the standard DNLS equation the 
background amplitudes have a monotone increasing behavior with u, see fig- 
ure (3. a). High frequency solutions cannot be continued up to that limit, 
and the continuation stops for values of v < 1. This result correlates well 
with the collective variable (particle perspective) predictions [19,25] where the 
non-persistence of travelling solutions is related to the growth of the Peierls- 
Nabarro barrier. In this respect, one observes a sudden increase in the back- 
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Fig. 2. Power Spectrum, (a) S(u) for the background of a 1/1 breather with 
Ldb = 5.057 at v = 0.05. In (b) formula (9) gives the contribution of a unique 
phonon (j = 0) with which agrees (eq. (8)) with results given by S(u>). (c) S(uj) for 
the background of a 1/2 breather with oj^ = 2.3842 at v = 1.00. In (d) formula (9) 
gives the contribution of seven phonons (j = 0, .., 6), but only five of them (0 — 4) 
are visible on S(u). Note that the amplitudes \Aj\ differ by orders of magnitude. 

ground amplitude near the continuation border. This result reinforces the 
interpretation of the background as an energy support to the core for sur- 
passing the (nonintegrable) Peierls-Nabarro barriers to mobility, and so its 
unavoidable presence for the existence of mobile breathers in the absence of 
integrability. 

The role of the background in the localized core mobility can be analysed as 
follows. As the solution is unambiguosly found to be $ = $ core + $ bcfc ^ the 
energy H, equation (5), of a mobile breather can be written as 



H = H[$ core ] + H{& ckg ] + H 



int 



(10) 



where H mt is the interaction energy, i. e. the crossed terms of $ core and Q hckg 
in the Hamiltonian. In the simplest case in which the background has a single 
resonant plane wave, its energy is a constant of motion (along with the total 
energy), so one obtains 

dH[¥ ore ] dH int 
dt ~ dT 
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(a) 



(b) 
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Fig. 3. (a) Background amplitude of different 1/1 resonant breathers as a function 
of the nonintegrability v. The amplitudes are zero in the A-L integrable limit (y = 0) 
and have a monotone increasing behavior with v. The value v = 1 corresponds to 
the standard DNLS equation, (b) Plot of H core of a 1/1 resonant breather with 
Ub = 5.056 as a function of the localization center xo for different values of v (0.04, 
0.08, 0.12, 0.16, 0.20, 0.24, 0.25 and 0.2512 (end of the continuation)). The amplitude 
of the oscillation of H core grows with v. (The minimum value of H core has been set 
to zero in order to compare the differents functions.) 



i. e. the variations of the core energy along the motion are balanced by the vari- 
ations in the core-background interaction energy. Equation (11) dictates the 
dynamics of any (eventual) effective (collective) variables intended to describe 
the mobile core in a particle-like description of the breather. 

One can compute the core energy variations directly from the numerical inte- 
gration of a solution by substracting (at each time step) the background from 
it. In figure (3.b) we plot the evolution of the core energy as a function of the 
core localization center, x$, denned by using the norm N (equation (6)) of the 
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Salerno model (2) as 

_ E re nln(l + /x|$r e | 2 ) ri9 x 
*o - ^ • (12) 

One observes that the core has extracted the maximum energy from the in- 
teraction term when the core passes over xq = n ± 1/2 (maxima of the PN 
barrier) and has given it back to the interaction term at xq = n (minima of 
the PN barrier). The bigger the PN barrier, the larger the interaction term 
(directly proportional to background amplitude) is. This result illustrate the 
role of the resonant background on the core mobility and the interpetation of 
its amplitude increase with v. The increase of nonintegrability, and the subse- 
quent growth of the PN barrier, demands aditional support of energy from the 
interaction term, which is achieved by an increase of the background amplitude. 



5 Conclusions 

We have used a (SVD)-regularized Newton algorithm to continue mobile dis- 
crete breathers in the Salerno model from the integrable Ablowitz-Ladik limit. 
Our results indicate that, away from integrability, a description of these so- 
lutions based exclusively on localized (collective) variables is incomplete. The 
solutions are composed by a localized core and a linear superposition of plane 
waves, the background, whose amplitudes differ orders of magnitude. The 
background plays an important role in the translational motion of the local- 
ized core. Exact mobile localization only exist over finely tuned extended states 
of the nonlinear lattice. Mobile "pure" (i.e. zero background) localization must 
be regarded as very exceptional (integrability). 
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